Current-voltage characteristics of tunnel Josephson junctions with a ferromagnetic interlayer 
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We present a quantitative study of the current-voltage characteristics (CVC) of diffusive superconductor/ 
insulator/ ferromagnet/ superconductor (SIFS) tunnel Josephson junctions. In order to obtain the CVC we 
calculate the density of states (DOS) in the F/S bilayer for arbitrary length of the ferromagnetic layer, using 
quasiclassical theory. For a ferromagnetic layer thickness larger than the characteristic penetration depth of 
the superconducting condensate into the F layer, we find an analytical expression which agrees with the DOS 
obtained from a self-consistent numerical method. We discuss general properties of the DOS and its dependence 
on the parameters of the ferromagnetic layer. In particular we focus our analysis on the DOS oscillations 
at the Fermi energy. Using the numerically obtained DOS we calculate the corresponding CVC and discuss 
their properties. Finally, we use CVC to calculate the macroscopic quantum tunneling (MQT) escape rate 
for the current biased SIFS junctions by taking into account the dissipative correction due to the quasiparticle 
tunneling. We show that the influence of the quasiparticle dissipation on the macroscopic quantum dynamics of 
SIFS junctions is small, which is an advantage of SIFS junctions for superconducting qubits applications. 

PACS numbers: 74.45.+C, 74.50.+r, 74.78.Fk, 75.30.Et 



I. INTRODUCTION 

The possibility to switch the ground-state of a Josephson 
junction from a to a tt phase state and the possible ap- 
plication of such junctions in quantum information led to a 
renewal interest in the study of the so called n Josephson 
junctions. The existence of such a transition was predicted 
more than thirty years ago,' however due to technological re- 
quirements only recently it was observed. The realization of 
n Josephson junctions was achieved in superconductor/ fer- 
romagnet/ superconductor (SFS) j unctions i2r2i Microscopi- 
cally, S/F hybrid structures are characterized by an unusual 
proximity effect, with a damped oscillatory behavior of the 
superconducting correlations in the F layer (for a review see 
Refs. 22-24 and references therein). This unusual proximity 
effect in S/F layered structures leads to a number of striking 
phenomena like the nonmonotonic dependence of their criti- 
cal temperature and the appearance of oscillations of critical 
current in SFS Josephson junctions as a function of the F layer 
thickness. ^'^ In particular the change of sign of the critical cur- 
rent corresponds to the so-called — ;r transition. 

On the other hand, SFS junctions, as any metallic junction 
exhibit very small resistances and therefore are not quite suit- 
able for those applications, for which active Josephson junc- 
tions are required. This problem can be solved by adding an 
additional insulating (I) layer to increase the resistance. SIFS 
junctions represent an interesting case for practical use of n 
Josephson junctions. For instance, a SIFS structure offers the 
freedom to tune the critical current density over a wide range 



and at the same time realize high values of the product of 
the junction critical current and its normal state resistance 
R„J^~2^ In addition, Nb based tunnel junctions are usually 
underdamped, which is desired for many applications. Due 
to these advantages, SIFS n junctions have been proposed as 
potential elements in superconducting classical and quantum 
logic circuits."^'"'' For instance, SIFS junctions can be used as 
complementary elements (tt - shifters) in RSFQ circuits (see 
Ref. 27 and references therein). Finally, SIFS structures have 
been proposed for the realization of so called ^-junctions with 
a (p drop in the ground state, where < (p < n?^ The prop- 
erties of SIFS junctions have been intensively studied both 
experimentally^— and theoretically.i^Srli However, proper- 
ties of the quasiparticle cuiTent have received relatively little 
attention so far, although they can be very important for the 
description of SIFS junctions as possible elements of super- 
conducting logic circuits. 

The purpose of this work is to provide a quantitative model 
describing the behavior of quasiparticle cuiTent in SIFS junc- 
tions as a function of parameters characterizing material prop- 
erties of the ferromagnetic interlayer. We also focus our study 
on the properties of the density of states (DOS) in S/F bi- 
layers and discuss the oscillations of the DOS at Fermi en- 
ergy. Finally, we calculate the macroscopic quantum tunnel- 
ing (MQT) escape rate for current-biased SIFS junctions by 
taking into account the dissipative correction due to the quasi- 
particle tunneling. Based on this we conclude that the influ- 
ence of the quasiparticle dissipation on the macroscopic quan- 
tum dynamics of SIFS junctions is small, which is an advan- 
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FIG. 1: (Color online) Geometry of the considered system. The 
thickness of the ferromagnetic interlayer is df. The transparency 
of the left (right) S/F interface is characterized by the coefficient 
YBXiBl) - The left interface is an insulating barrier, y^l ^ 1 (shown by 
a black line), while the right interface is transparent, Jbi <S 1 (shown 
by a grey line). 

tage of SIFS junctions for quantum logic (qubit) applications. 

The paper is organized as follows. In the next section we 
formulate the theoretical model and the basic equations. In 
Sec. |III]we solve the nonlinear Usadel equations numerically 
for arbitrary length of the ferromagnetic layer, and calculate 
the DOS in the F layer. We compare these results with an an- 
alytical expression for DOS in case of a long SIFS junction, 
i.e. when the thickness df of the ferromagnetic layer is much 
larger than the decay length of the characteristic supercon- 
ducting correlations in the ferromagnet ^ji. We also discuss 
the oscillations of the DOS at the Fermi energy. In Sec. |IV] 
we present the current-voltage characteristics of SIFS junc- 
tions for different parameters of the ferromagnetic interlayer. 
In Sec.|V] we use these data to calculate the MQT escape rate 
for current-biased SIFS junctions by taking into account the 
dissipative effect of the quasiparticle tunneling. Finally we 
summarize the results in Sec. |VI] 



II. MODEL AND BASIC EQUATIONS 

We consider a SIFS junction such as the one depicted in 
Fig. [T] It consists of a ferromagnetic layer of thickness df 
and two thick superconducting electrodes along the x direc- 
tion. The left and right superconductor/ ferromagnet inter- 
faces are characterized by the dimensionless parameters Jbi 
and Yb2?^'^^ respectively, where 7^1,^2 = ^Bi,B20-n/^«, Rb\,b2 
are the resistances of the left and right S/F interfaces, respec- 
tively. On is the conductivity of the F layer, ^„ = ^jDf /2nTc, 
Df is the diffusion coefficient in the ferromagnetic metal and 
Tc is the critical temperature of the superconductor (we as- 
sume h — kg — 1, except for Sec. |V]l. We also assume that 
the S/F interfaces are not magnetically active. We will con- 
sider the diffusive limit, in which the elastic scattering length 
£ is much smaller than the decay characteristic length ^fi = 
min [^^■i|,<^^i|] (for the definitions of lid) see Eqs. (fTSIl be- 
low). 

We assume that the tunneling barrier is located at the left 
S/F interface, while the right interface is perfectly transparent, 
this means that Jbi ^ 1 while 7^2 <C 1. In this case the left S 
layer and the right F/S bilayer in Fig. [T] are decoupled and we 
can calculate the quasiparticle current through a SIFS junction 
using the standard tunneUng formula^ 

I^^J dENsiE-eV)NfiE)[f{E-eV)-f{E)], (1) 



where Ns{E) = \E\@{\E\ - A) /y^E^ - is the BCS DOS, 
&{x) is the Heaviside step function, Nf{E) is the DOS in 
the ferromagnetic interlayer at x = —df/2, f{E) = [1 + 
exp(£'/r)]^' is the Fermi function, and/? = /?Bi. BothNs{E) 
and Nf{E) are normalized to their values in the normal state. 
In particular at zero temperature, T = 0, the current acquires 
the form, 

I reV-A 

I = @{eV - A) — / dEN, {E - eV)Nf{E) . (2) 
eR Jo 

To obtain Nf{E) we notice that since Jbi ^ 1> the left su- 
perconducting lead does not influence the DOS in the ferro- 
magnetic interlayer (to zero order in the barrier transparency). 
This reduces the problem to the following: we need to find the 
DOS of a single F/S bilayer, which can be done by solving the 
Usadel equations in the ferromagnetic layer. 

Using the (J-parameterizations of the normal and anoma- 
lous Green functions, G = cos 0, F = sin 0, we can write the 
Usadel equations in the F layer as^"— 

Dfd^ef^(n ( 1 \ 

2 ^\^^ + ^ ) 

+ 1 sin (0^1 + 0^;) ± -1 sin (0^t - %) . (3) 

where the positive and negative signs correspond to the spin 
up t and spin down \, states respectively. In this notation 
the spin up state corresponds to the anomalous Green func- 
tion F| ^ {xj/^Xj/i) while the spin down state corresponds to 
F| ^ where V't(i) '^^^ electron fermionic opera- 

tors. The CO — 2nT{n + j) are the Matsubara frequencies, 
and h is the exchange field in the ferromagnet. The scatter- 
ing times are labelled here as Tj, Xx and Tjo, where T,(^) cor- 
responds to the magnetic scattering parallel (perpendicular) 
to the quantization axis, and 1^0 is the spin-orbit scattering 
time.i^-^ 

We consider here ferromagnets with a strong uniaxial 
anisotropy, in which case the magnetic scattering does not 
couple the spin up and spin down electron populations, i.e. 
the perpendicular fluctuations of the exchange field are sup- 
pressed (T^r' ~ 0). Therefore, we will neglect ix in our con- 
sideration and denote as a magnetic scattering time T„,. We 
will also consider ferromagnets with weak spin-orbit interac- 
tions and henceforth also neglect the spin-orbit scattering time 
Tso- In this case the Usadel equations in the ferromagnetic 
layer for different spin projections are not coupled any more 
and can be written as, 

Df ^'Q/td) _ f^,^ , , £^!%m)^ sin (A\ 

T^^" (,''^''' + ^;;~j'^"^/ta)' (4) 

while in the S layer the Usadel equations take the form 
D d^O 

— ^ -T-^ = (0 sin 0s — A(x) cos 0j . (5) 
2 dx'^ 

Here is the diffusion coefficient in the superconductor and 
A{x) is the superconducting pair potential. Notice that in the 
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latter equation we have omitted the subscripts 't (i)' because 
both equations are identical in the superconductor 

Eqs. (|4|, (|5]) should be complemented by the self- 
consistency equation for the superconducting order parameter 
A, 

A(^) In I = Trr ^ - sin 0,^ - sin , (6) 

and by the boundary conditions at the outer boundary of the 
ferromagnet 



dx 



and at the F/S interface,— 



-df/i 



clfl2 



0, 



dd, 



(7) 



^"^( dx ) ,i..n V dx 



39 



^«7b2(^) - sin (0,-0/) ^,^/2' 

df/2 



df/2 
df/2 



(8a) 
(8b) 



where 7= ^i C7„/4,<7j , C7, is the conductivity of the S layer and 
£,s = ^Ds/ItzTc. The parameter 7 determines the strength 
of suppression of superconductivity in the right S lead near 
the interface compared to the bulk: no suppression occurs for 
7=0, while strong suppression takes place for 7 ^ 1 . In our 
numerical calculations we will assume small 7 <C 1 . Notice 
that the interface parameters do not depend on the spin di- 
rection. In other words we are not considering spin-active 
interfaces. In the case of spin-active barriers, one should 
use the boundary conditions introduced in Refs. 42-44, rather 
than the standard Kupriyanov-Lukichev boundary conditions 
[Eqs. ©]. 

To complete the boundary problem we also set a boundary 
condition atx = °o 



0,(00) = arctan 



CO 



(9) 



where the Green functions acquire the well known bulk BCS 
form. Eqs. dUi-® represent a closed set of equations that 
should be solved self-consistently. As it will be discussed in 
the next section, the knowdlege of the Green function will al- 
low us to compute the DOS at the outer F boundary. 



III. DENSITY OF STATES IN THE F/S BILAYER 

The DOS Nf{E) normalized to the DOS in the normal state, 
can be written as 



Nf{E)^[Nf^{E)+Nf^{E)] /2, 



(10) 



where Nf^(^^^^{E) are the spin resolved DOS written in terms 
of spectral angle 0, 

%a)(^) =Re[cos0/^(;)(/a)^£ + /O)] . (11) 



To obtain Nf, we use a self-consistent two-step iterative 
procedure^^'^-"^. In the first step we calculate the pair po- 
tential coordinate dependence A(.ji:) using the self-consistency 
equation in the S layer, Eq. (|6]l. Then, by proceeding to the 
analytical continuation in Eqs. (|4]i, (|5]l over the quasiparticle 
energy ico ^ E + iO and using the A(jic) dependence obtained 
in the previous step, we find the Green functions by repeating 
the iterations until convergency is reached. 

Before showing the numerical results we consider an ana- 
lytic limiting case. If the F layer is thick enough (df ^ ^/i) 
and 7 = in Eq. dSal l. the DOS at the free boundary of the 



ferromagnet can be written ai 



,29.48 



/t(i) 



(£)=Re[cos0fo|(;)]«l--Re0j 



Here 06|(|) is the value of 9f at x - 



2"- b^il)- 

df/2, given by 







8F(£) 



v/(l-7?2)^'2(£) + l + l 



exp 



df 



(12) 



(13) 



where £,f = yjDjjh. In Eq. ([13]) we use the following nota- 
tions, 

^/2Jh^J-iER±ih^\|x^, (14a) 
(l/T,„)H£ff±//2 + l/T„,)"', (14b) 
A 



F{E) = 



-ier + Ja'^ 



Er = E + iO. 



-£■2 



(14c) 



Here, we again adopt the convention that a positive (negative) 
sign in front of h corresponds to the spin up state '[ (spin down 
state 1). Hereafter we will write spin labels t (|) explicitly 
only when needed. 

From Eqs. (|l2li-(|T3]l we obtain for the full DOS the follow- 
ing expression in the limit df ^ ^fi. 



l6FHE)exp[^p'f 

t^i(V(i-rj2)F2(£)TT+i)2 



(15) 



At this point, we define the characteristic decay and oscillation 
lengths ^/i.2ta) as 



'/ltd) 



'/2t(i) 



1 



'E^h 



1 



'E^h 



1 



(16a) 



(16b) 



(16c) 



In the absence of m agnetic sc attering both lengths coincide 
and are equal to £,fy^h/\E^h\ for different spin orientations. 
One can rewrite Eq. ([TsT i in the following form. 



^exp 

t4 



Scos X 



2df 

'in 

2df 

^W2 



Asin X 



2df_ 

if2 



(17) 
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FIG. 2: (Color online) DOS Nf{E) on the free boundary of the 
F layer in the F/S bilayer calculated numerically in the absence of 
magnetic scattering (1/t„A = 0) for different values of the F layer 
thickness df, h/h = 4,T = 0.17^. Parameters of the F/S interface are 
y= Ybi = 0.01. (a): df/^„ = 0.5, (b): df/^„ = 1, (c): d^^„ = 2, (d): 



df/^„ = 3. The approximate analytical solution, Eq. l llSb . is shown 
by dashed red lines. 

where the coefficients A, 23 and X can be obtained by ex- 
pansion of the real part in Eq. ( fTSl ); only two of them are 
independent. This form explicitly shows the damped oscil- 
latory behavior of superconducting correlations in the F layer 
The lengths ^fi,2 also the lengths of decay and oscilla- 
tions of the critical current in SIFS junctions (see Eqs. (26) in 
Ref.l29l). The period of the DOS oscillations is approximately 
twice smaller than the period of the critical current oscillations 
and the exponential decay is approximately twice faster than 
the decay of the critical current»2E 

Now we turn to the exact numerical solution. The obtained 
energy dependencies of the DOS at the free F boundary of the 
F/S bilayer are presented in Figs.|2][3]|4] The exchange field is 
chosen such that h> A, which corresponds to the experimental 
situation. 

Figure|2]shows the DOS energy dependence for different df 
in the absence of magnetic scattering. At small df we observe 
the DOS double peak due to the Zeeman splitting of the BCS 
peak at £■ = A. Most probably in the experiments, the BCS 
Zeeman splitted peak as presented in Fig. |2](a) will be seen 
as a single peak due to many-body interaction effects, which 
introduce a finite lifetime (damping) of the quasiparticles. We 
also observe that at small df and relatively small exchange 
field h, full DOS turns to zero inside a minigap, which van- 
ishes with the increase of df. 

The minigap also exists in the normal metal (N) DOS in 
the S/N bilayers. If the thickness d„ of the normal metal is 
larger than the coherence length, the characteristic scale of the 
minigap is set by the Thouless energy, Eji, — D„/d^, where 
D„ and is the diffusion coefficient of the normal metal. In the 
F layer of the S/F bilayer, the exchange field li shifts the DOS 
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FIG. 3: (Color online) DOS Nf{E) on the free boundary of the 
F layer in the F/S bilayer calculated numerically in the absence of 
magnetic scattering (1/t,„A = 0) for different values of the exchange 
field /j. Parameters of the F/S interface are y= 7^2 = 0.01, r = 0.irc. 
Plots (a) and (b): df/£,„ = 1; plots (c) and (d): df/£,„ = 3. For plots 
(a) and (c) solid black line corresponds to h /A = 2, dashed red line 
to h/A = 2.5, dash-dotted blue line to h/A = 3. For plots (b) and (d) 
solid black line corresponds to /7/A = 4, dashed red line to /j/A = 5, 
dash-dotted blue line to /j/A = 6. 

for the two spin subbands in opposite directions, therefore the 
critical value he of the exchange field at which the minigap in 
the spectrum closes can be roughly estimated 



he 



ETh^Df/dj. 



(18) 



This equation shows the qualitative tendency: for smaller df 
a higher h is needed to close the minigap [see also Fig.[3](a)]. 
The estimation, Eq. ( fTSI l, is only valid in the absence of mag- 
netic scattering, since t,„ also influences the minigap^ [see 
alsoFig.|4](a)]. 

In Fig. |2] we also observe that after the minigap closes, 
the DOS at the Fermi energy Nf{Q) rapidly increases to val- 
ues larger than unity with further increase of df; then it os- 
cillates around unity while its absolute value exponentially 
approaches unity [see also Fig. |5l. This is the well-known 
damped oscillatory behavior of the DOS in F/S bilayers. Ex- 
perimental evidence for such behavior was provided by Kon- 
tos et a/i— In the case of long enough ferromagnetic layer we 
also observe the DOS peak at E — h, which was previously 
discussed in Ref. '52'. A similar effect was also discussed in 
N/F/S structures, where it was shown that a zero energy peak 
appears in DOS if Eji, ~ 

We also show in Fig. |2] the analytical approximation, 
Eq. ( fTSl ), which is in good agreement with the numerical re- 
sult for thick enough feiTomagnetic layers. In the numerically 
obtained curves the peak atE = his smeared because of finite 
7 = 0.01 for the transparent F/S interface atx = df/2. 

In the absence of magnetic scattering we rewrite the ana- 
lytical DOS expression, Eq. ( fTSl ), for £ > A in the following 
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FIG. 4: (Color online) DOS Nf{E) at the free boundary of the F 
layer in the F/S bilayer calculated numerically for = 1/ TmA = 0.5 
(solid black line), a,„ = 1 (dashed red line), and a,„ = 3 (dash-dotted 
blue line) for different values of the F layer thickness rfy, h = 4A, 
T = O.irc- Parameters of the F/S interface are y = 7^2 = 0.01. (a): 
dflt,n = 0.5, (b): dfj^n = 1, (c): dfj^n = 2, (d): rf//<^„ = 3. For 
plots (c) and (d) the curves with = 3 are not shown since they 
are of the order of unity at corresponding scale. Inset of the plot (a): 
Nf{E) dependence for df/^„ = 0.5 for higher values of a,„; a,n ~ 5 
(solid black line), a,„ = 7 (dashed red line), a,„ = 15 (dash-dotted 
blue line). 



way, 



16A2cos ( ^ 



Nf{E) = 1 - 



\E±h\ 



\E±h\ 



-e)(v^ 



2e)2 



(19) 



where e = \/E'^ — A^. We can clearly see the exponential 
asymptotic of the peak atE = h from the Eq. (fT9] l. We should 
keep in mind that Eq. ( fT9] l is valid for large df/^f, but never- 
theless we may qualitatively understand why we do not see the 
peak at E — h for small ratio of df/^j: if this factor is small 
the variation of the exponent e\p{~2{df/^f)y^\E — h\/h} 
near the point E ~ h is also small. The peak is observable 
only for h of the order of a few A. For larger exchange fields 
the peak is very difficult to observe, since the energy depen- 
dent prefactor of the exponent in Eq. ( fT9] ) decays as E^^ for 
£>A. 

Figure [3] shows the DOS energy dependence for different 
values of the exchange field h in the absence of magnetic 
scattering. For stronger exchange field the minigap closes at 
smaller df, in qualitative correspondence with Eq. ( fTSl l. From 
numerical calculations we obtain the following condition [see 
also Fig.|5](c)], valid for T,;;' ^ 0,^- 
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FIG. 5: (Color online) Dependence of SNf as a function of the F 
layer thickness df for different exchange fields (a) and magnetic scat- 
tering times (b); dependence of SNf as a function of exchange filed 
(c) and magnetic scattering time (d) for different df. The tempera- 
ture r = O.lTc. Parameters of the F/S interface are y = 7^2 = 0.01. 

(a) : no magnetic scattering (1/t,„A = 0), h/A = 2 (black solid 
line), h/A = 4 (red dashed line), h/A = 6 (blue dash-dotted line); 

(b) : h/A = 4, a,,, = 1/t:„,A = (black solid line), a,„ = 1 (red 
dashed line), «„, = 3 (blue dash-dotted line); (c) no magnetic scat- 
tering, df/^n = 0.5 (black solid line), df/^,, = 1 (red dashed line), 
df/£,n =2 (blue dash-dotted line), df/^n = 3 (black short-dashed 
fine); (d) h/A = 4, df/£,„ = 0.5 (black solid line), df/^„ = 1 (red 
dashed line), df/£,„ = 2 (blue dash-dotted line), df/^n = 3 (black 
short-dashed line). 

In Fig. Owe also observe the peak atE = h; at large enough 
exchange fields its amplitude can be much larger than the am- 
plitude of the peak at £ = A (see Fig.[3](d), blue dash-dotted 
curve). The existence of the DOS peak at E = h gives a pos- 
sibility to measure the exchange field directly in experiment 
by measuring the F/S bilayer DOS in compounds with small 
magnetic scattering (since magnetic scattering is smearing the 
peak, see below). For example, in Ref. were reported ex- 
change fields for Pdi_;fNiT with different Ni concentration, 
obtained by a fitting procedure. Considering Nb as a super- 
conductor with A =1.3 meV, we can estimate the exchange 
field in Pdi-jcNi^: for 7% of Ni fitting gives h =2.8 meV, 
which is 2.2A, and for 11.5% of Ni h =3.9 meV, which is 
3A.^^ It is interesting to use direct measurements of the DOS 
peak at E = I1 to check these fitting predictions of Ref. [Hi 

Ferromagnetic metals with exchange fields of the order of 
few A are crucially important for the fabrication of SIFS junc- 
tions, valid for superconducting logic applications. Presently 
used ferromagnets have h ^ A, and therefore short oscilla- 
tion length [see Eq. (I16cl )l. which makes it difficult to control 
the F layer thickness. In already existing SIFS structures the 
roughness is often larger than desired precision of df.^^. We 
hope that our results will trigger the experimental activity in 
finding ferromagnetic alloys with /; of the order of few A. 

Figure |4] shows the DOS energy dependence for different 
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values of magnetic scattering time. Similarly to Fig. [3] for 
stronger magnetic scattering the minigap closes at smaller df. 
Also the DOS peak at £ = h, visible for long enough fer- 
romagnetic layer, is smeared. The analytical solution (not 
shown), Eq. ( fTSl ) also agrees quite well with the numerical 
results for t// ^ ^/i . 

Although our results are obtained for weak ferromagnets, 
they can in certain cases be extended for ferromagnets with 
strong exchange fields, h ^ A. In the absence of mag- 
netic scattering the Usadel equation in energy representation, 
Eq. (HJi, can be rewritten as. 



fm 



2 dy^ 



h 



=F 1 sin 



/t(i)' 



(21) 



where y — x/ is the dimensionless coordinate. In the case 
of /; ^ A, we can neglect the first term on the right hand side 
of Eq. ( |2TI ) to obtain the subgap DOS. Thus, in that case the 
subgap structure scales with the length and for example the 
results presented in Fig.|2]for h — AA also describe the DOS in 
the case of a high exchange field if one scales df correspond- 
ingly. This procedure does not apply however for h = 2A (see 
Fig.O, since in that case one cannot simply neglect the term 
£//iinEq. (EB- 

To show explicitly the aforementioned DOS oscillations at 
the Fermi energy,^' in Fig.|5]we plot the numerically calcu- 
lated function 

5Nf{dfX^,n) = \\-NfQ\, Nfo=Nf{E = 0). (22) 

Using Eqs. (fT2l l-(fT7]l and (l22T l we get the analytical expres- 
sion for the function dN, valid for df '3> ^fi. 



5N = 32 



Re 



(./2-77o' + l)2 



exp -po 



where 



Po = \/2/hy^ ih + I / 

J70 = (l/Tm)(//i+l/T,„)"^ 

At vanishing magnetic scattering we obtain. 



dN: 



32 



3+2V2 



2df 

cos I ) exp 



2df 

17 



(23) 



(24a) 
(24b) 



(25) 



in which case the characteristic lengths of decay and oscilla- 
tions are equal to ^f. 

The dependence of 5Nf{df) on the ferromagnetic layer 
thickness df at different values of exchange field and mag- 
netic scattering time is presented in Fig.|5](a) and (b). From 
Fig.|5](a) we can see that with increasing exchange field h the 
minigap closes at smaller df in agreement with Eq. ( |20| ). the 
period of the DOS oscillations at the Fermi energy decreases, 
and the damped exponential decay occurs faster. This is easy 
to see from Eq. ( |25] |. since in the absence of magnetic scatter- 
ing 8N depends on h only as a function of 

From Fig.|5](b) we can see that with increasing a,„~l/ T„,A 
the period of the DOS oscillations on the contrary increases. 



although the minigap also closes at smaller df and the damped 
exponential decay occurs faster To understand this behavior 
we rewrite here the decay and oscillation lengths, Eq. ( fTSI l. at 
the Fermi energy. 



1 

1 




(26a) 



(26b) 



We see from these equations that with increasing a,„ the 
length of decay ^fi decreases, while the length of oscillations 
^f2 increases. 

The dependencies of 5Nf on exchange field and magnetic 
scattering time are presented in Fig.|5](c) and (d), correspond- 
ingly. In Fig.|5](c) we see oscillations of the DOS at the Fermi 
energy 8Nf{h) around unity with increasing exchange field in 
the absence of magnetic scattering. In Fig.|5](d) we show the 
function 5A^(t„,). It is interesting to note that its behavior can 
be both oscillatory and also monotonous. When the parame- 
ter a,„ increases starting from the minigap state (black solid 
curve) it is totally monotonous: increasing a,„ the minigap 
closes and the DOS starts to increase to unity, but never over- 
shoots unity [we checked this up to a,„ ~ 80, see also Fig. |4] 
(a)]. If we start from the state where the minigap is already 
closed, we first observe oscillations, but then again a switch 
to monotonous behavior For intermediate F layer thicknesses 
we see just one oscillation and then DOS monotonously ap- 
proaches unity (we checked this up to a,„ — 80), while for 
thicker feiTomagnets (df/^n = 3) we observe two oscillations 
and then a monotonic behavior. 

The dependencies 5N{h) and 5N{t,„) can be important if 
in the experiment the material properties of the ferromagnetic 
interlay er, i.e. exchange field h and magnetic scattering time 
Tm can vary with some external parameter, for example tem- 
perature, magnetic field, etc. 

Before turning to the calculation of the CVC we discuss 
briefly a recent experiment^^ in which a pronounced double 
peak in the DOS of Ni/Nb bilayers were reported. This dou- 
ble peak cannot be explained within our model based on the 
Zeeman splitting. The reason for the double peak in Ref. fs^ 
remains controversial. In Ref. 113 it was numerically fitted by 
adding an extra parameter to the model, characterizing spin- 
active interfaces. However, this fit is far from being satisfac- 
tory. Nevertheless there is another feature of the DOS ob- 
served in Ref. which can be explained within our model: 
by increasing df the "normal" peak at £ = A [which is the 
BCS Zeeman split peak in Fig. |2la)] is "inverted" [Fig. |2lb) 
and (c)] and becomes "normal" again [Fig.|2td)] as df is fur- 
ther increased. According to our model, at £ = A Eq. ( fT9] l 
reduces to the following expression, 



' 2d 

Nf{A) = l + 16^cos( -j^ 




(27) 



This expression explains the inversion of the peak at £ = A 
as a function of df. The peak is "normal" ("inverted") if the 
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FIG. 6: (Color online) Current- voltage characteristics of a SIFS junc- 
tion in the absence of magnetic scattering for different values of the 
F-layer thickness d j . The temperature T = Q.\Tc. The exchange field 
/! = (black line, which corresponds to the case of a SINS junction), 
h/h = 2 (blue dash-dotted line), and /i/A = 4 (red short-dashed line), 
(a): df/^„ = 0.5; (b): df/^„ = 1, (c): df/t,,, = 2, and (d): df/t,„ = 3. 
Insets in (a), (b) and (c) are explained in the text. 

DOS at A is larger (smaller) than unity. This variation is due 
to the sign of the cosine function in Eq. ( |27] ). which depends 
on the df/^f ratio. 



rV. CURRENT- VOLTAGE CHARACTERISTICS OF A SIFS 
JUNCTION 

In this section we calculate the current-voltage characteris- 
tics (CVC) of a SIFS junction at low temperature, T — O.IT^, 
using Eq. ([T]i and DOS Nf{E,df,h,T,„) numerically obtained 
in the previous section. 

Figure |6] shows the CVC of a SIFS junction in the absence 
of magnetic scattering. For comparison we also present the 
CVC of a SINS tunnel junction, i.e. a junction with a normal 
metal interlayer instead of a ferromagnet (h = 0). SINS struc- 
tures were studied previously in Ref. 46. We observe several 
features of SIFS CVC which are the signatures of the proxim- 
ity effect in the S/F bilayer 

For thin enough F layer we observe the "kink" on the CVC 
at eV « 2A [Fig. |6](a)], which corresponds to the case when 
the DOS Nf{E) exhibits a pronounced minigap. The corre- 
sponding DOS energy dependence (/i/A = 4, df/^„ = 0.5) is 
shown in the inset. We can also see that for a certain range 
of parameters the CVC of a SIFS junction exhibit a non- 
monotonic "wave" behavior We can observe it for h/A = 4- 
(red short-dashed line) in Fig. |6](b) and for h/A = 2 (blue 
dash-dotted line) in Fig. |6](c). This behavior corresponds to 
the case when the DOS Nf{E) minigap is already closed and 
the Nf{Q) at the Fermi energy is larger than unity. The cor- 
responding DOS energy dependencies are presented in the in- 
sets of the plots |6](b) and (c). 




FIG. 7: (Color online) Current- voltage characteristics of a SIFS junc- 
tion for l/tmA = 1 (black solid lines) and 1/t,„A = 3 (red short- 
dashed lines) for different values of the F-layer thickness dj. The 
temperature T = O.lTc. The exchange field h/A = 2: plots (a) and 
(b) and h/A = 4: plots (c) and (d). The thickness df/£,„ = 0.5: plots 
(a) and (c) and dj/^„ = 1: plots (b) and (d). Insets in (a) and (c) are 
explained in the text. 

At large enough df and exchange fields the DOS Nf{E) w 1 
and the current, Eq. ([T]l, is given by the same equation as the 
current in the NIS tunnel junction. 



1 



dEN,{E)[f{E-eV)~fiE)]. 



(28) 



At small temperature T <C 71-, this equation is well approxi- 
mated by taking T — 0, 



1 



I = @ieV^A) — 

eR J A 



eV 



EdE 

V£2-A2 



@(eV-A) 



v/(ey)2-A2 



(29) 



The red short-dashed line in Fig. |6] (d) {df/^„ = 3, h/ A = 
4, no magnetic scattering) almost coincides with this result 
except for the small region eV ~ A, since for our numerically 
calculated curves we fix temperature T = O.lTc <C T^. 

Figure [T] shows current-voltage characteristics of a SIFS 
junction in case of finite magnetic scattering. Here for thin 
F layers we observe a "double-kink" structure, see Fig.|2](a) 
and (c). It corresponds to the DOS Nf{E) with small mini- 
gap and finite subgap value smaller than unity. Such a DOS 
structure is typical in the presence of magnetic scattering and 
thin enough ferromagnetic interlayer, see Fig.|4](a). The cor- 
responding DOS energy dependencies are presented in the in- 
sets of the plots |6](a) [1/t„,A 3,/z/A = 2,<i//4, = 0.5] and 
(b) [1/t„,A 3, h/A = A,df/^„ ~ 0.5]. For finite magnetic 
scattering the non-monotonic features of CVC are smeared. 
We do not show the curves for df/^„ > 2, since they do not 
significantly differ from the curves obtained from Eq. ( |29] ). 

Figures|6](a) and|2](a) show that for thin enough ferromag- 
netic layer the current has an onset in the interval [A, 2 A] (for 
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fective action 5eff is given by 



FIG. 8: (Color online) (a) Schematic of a current biased SIFS Joseph- 
son junction, /ext is the external bias current, (b) Potential U{^) 
v.s. the phase difference between two superconductors. cOp is the 
Josephson plasma frequency of the junction. 

temperatures T <^ Tc). The value of this onset, according to 
Eq. (|2]i, is A + Eg, where Eg is the DOS minigap, < Eg < A. 
Increasing exchange field, magnetic scattering and/or F layer 
thickness, the minigap closes and the current turns to zero 
at eV < A, having an onset at eV — A. The dependence of 
the minigap Eg on the parameters characterizing the material 
properties of the ferromagnetic interlayer is discussed in Sec- 
tionUn] 

We conclude that we observe interesting features in the 
SIFS CVC if the DOS Nf{E) near the insulating barrier has a 
nontrivial shape in the subgap region. In case when Nf ^ I, 
these features disappear and the CVC coincide with those of 
the NIS tunnel junction, Eq. 



MACROSCOPIC QUANTUM TUNNELING IN A SIFS 
JUNCTION 



In this section, motivated by experimental studies on the 
(thermal and quantum) switching^^ and quantum coherent os- 
cillations^^' in SFS and SIFS junctions, we calculate the MQT 
escape rate in a current biased SIFS junction as shown in 
Fig. [8] (a). The CVC obtained in the previous section enable 
us to investigate the influence of the quasiparticle dissipation 
on MQT 

It is important to note that MQT can be used as a mea- 
surement process of a superconducting phase qubiti^ Thus, 
the calculation of the MQT rate by taking into account the 
quasiparticle dissipation will be very important for analyzing 
the fidelity of measurement process for phase qubits. In the 
following calculation, for simplicity, we have ignored the in- 
fluence from an environmental circuit on MQT which can be 
experimentally reduced by a noise filtering technique.^- 

The partition function of a junction can be described by 
an imaginary-time functional integral over the macroscopic 
variable (the phase difference between two superconduc- 
tors),"-^'* i.e.. 



Z= / D0(T)exp 



(30) 



Ses[<P] 



np 



dx 



rhfl t-np 

Sa[(l)]^- d-i dx'a{T-t') 
Jo Jo 



M fd(j){T:) 



2 \ dx 



-U{(j)) +Sam, (31a) 
2 

(31b) 



In this equation, j3 = l/ksT, M = C{h./2ef' is the mass (C is 
the capacitance of the junction) and the potential t/(0) can be 
described by a tilted washboard potential [Fig.[8](b)], i.e.. 



U{(p) = -E J [sgn{Ic) cos (p+y^], 



(32) 



with y = Ie.xt/\Ic\, where Ej = h\Ic\/2e is the Josephson cou- 
pling energy, Ic is the Josephson critical current, and I^xt is 
the external bias current, respectively. The dissipation kernel 
a(T) is related to the quasiparticle current / under constant 
bias voltage V by 



' d(0 
e Jo 2n 



, 61-63 



exp(-a)T)/ V 



hco 



(33) 



at zero temperatureii 

As clearly seen from Figs. 6 and 7, the CVC has a gap 
structure due to the isotropic superconducting gap in left su- 
perconductor electrode. In such a case, the dissipation ker- 
nel a decays exponentially as a function of imaginary time 
T for |t| ^ Ti/A. The typical dynamical time scale of the 
macroscopic variable (j) is of the order of the inverse Joseph- 
son plasma frequency cOp = yj2e\Ic\/h.C{l — y^)'/^ which is 
much smaller than A. Thus, the phase varies slowly with the 
time scale given by h/A, and we can expand 0(t) — in 
Eq. ( I31bl i about T = t'. This gives 



5c r'^P 

2 Jo 



1 2 



2e dx 



where 



dC = 2 



dxa{x)x'^ 



(34) 



(35) 



Hence, the dissipation action Sa acts as a kinetic term so that 
the effect of the quasiparticles results in an increase of the 
capacitance, C — C + 5C = Cren- 

In case of a thin ferromagnetic layer (df ~ 0.5^„) we nu- 
merically obtain 5C ^h/AR for CVC presented in Figs.|6](a), 
|7](a), (c) ["kink" and "double-kink" structures]. For thick fer- 
romagnetic layer we can use Eq. ( |29] | to calculate 5C, 



4h 
nAR Jo 



dxx 



dz e' 



1 



2h 

AR' 



(36) 



For intermediate df we numerically find 5C ~ (1-^2) h/AR. 
Considering Nb as a superconductor (A = 1.3 meV) we there- 
fore obtain 5C ~ (0.5 pF, where r is the junction re- 
sistance R in £1. To insure a small dissipative correction of 
capacitance, 5C <C C we have a constraint. 



In the strong insulating barrier limit, i.e., Jex ^ 1, the ef- 



RC:^h/A, 



(37) 
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Without dissipation Tq 




With dissipation T 
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i? = 2mn 



^extj Pel 

FIG. 9: (Color online) The MQT escape rate for a current-biased 
SIFS junction as a function of the bias current /„f for several values 
of the junction resistance R. To (black solid line) and F (red dashed, 
yellow dot-dashed, and blue dotted lines) are the MQT escape rate 
without and with quasiparticle dissipation, respectively. Parameters 
are C = 800 pF, A = 1.3 meV, and = 500 ^iA.'^ 

i.e. the typical time constant RC of a SIFS junction should 
be much larger than the dynamical damping scale for the dis- 
sipation kernel (x{x). For example, in Ref. 13 the following 
parameters of a Nb/Al203/Nio.6Cuo.4/Nb SIFS junction were 
reported, C =800 pF and R =55 m£2, which correspond to 
5C w 10-:- 18 p F. Thus, even for the low resistive tunnel bar- 
rier in Ref. 13, we have the condition Eq. (|37] | satisfied. 

In order to see the effect of the quasiparticle dissipation on 
macroscopic quantum dynamics, we will investigate MQT in 
current-biased SIFS junctions. The MQT escape rate F from 
the metastable potential at zero temperature is given by^^ 



F = lim -- ImlnZ. 



(38) 



By using the Caldeira and Leggett theory,^** the MQT rate is 
approximated as 



F = ^Vl20;rB exp(-B), 
2% 



where 



(39) 



(40) 



is the renormalized Josephson plasma frequency and B — 
•Seffi^fi]/^ is the bounce exponent, that is the value of the ac- 
tion 5eff evaluated along the bounce trajectory ^b(t). The an- 
alytic expression for the bounce exponent is given by 



B 



12 
5^ 



le 



|/c|Qe„ (1 



(41) 



At high temperatures, the thermally activated decay domi- 
nates the escape process. Then the escape rate is given by the 
Kramers formula, F = (a)p/2;r)exp(— f/o/^B7")> where f/o is 
the barrier height.^^ Below the crossover temperature T* , the 



escape process is dominated by MQT. In the low dissipative 
(underdamping) cases, T* is approximately given by^^^ 



Inkn ■ 



(42) 



Here (y) = /q' dyP(y)y is the average switching current, where 
Pij) is the switching current distribution which is related to 
the escape rate F as^ 



PW = -F(y)exp 



1 



Y{y')dy' 



(43) 



In this equation, v = \dyldt \ is the sweep rate of the external 
bias current. Importantly, T* is reduced in the presence of 
dissipative effects. 

By using Eqs. ( [36] l and ( [39] l, we calculate F and compare it 
with the case without the quasiparticle dissipation. In Fig. 9, 
we numerically plot F and Fq for C = 800 pF and \lc\ = 500 
/xAfi^ and several values of R, where Fq is the MQT escape 
rate in the absence of the quasiparticle dissipation [Cren C 
in Eq. ([39]ll. As seen in this figure, F shows strong depen- 
dence on the junction resistance R, and F is almost identical 
to Fq in the case of large R, e.g., R = 55 mQ. (> h/CA w 0.63 
mi2) which corresponds to the actual SIFS junction.'^ We 
also calculate T* for a realistic case {R = 55 m£2) and found 
that T* = 7.4mK for the dissipative case (Cren = C + 8C) and 
T* = 7.5mK for the dissipation-less case (Cren — C)- As ex- 
pected, the T* suppression is small enough to allow experi- 
mental observations of MQT. Thus we can conclude that the 
influence of the quasiparticle dissipation on the macroscopic 
quantum dynamics of SIFS junction is very small for the case 
when the condition Eq. ( |37] ) is hold. This fact strongly sug- 
gests the great advantage of realistic SIFS junctions for qubit 
applications. The smallness of the quasiparticle dissipation in 
SIFS junctions is due to the superconducting gap in the left S 
electrode and the strong insulating barrier (7^1 ^ 1) between 
the left S and F layers. 

It is important to note that such a weak quasiparticle- 
dissipation nature of MQT has been also predicted in 
TT-junctions based on a S/ferromagnetic insulator (FI)/S 
junctions<^^j^ However no ferromagnetic insulator based 
Josephson junctions have been experimentally realized at 
present. On the other hands, the fabrication of SIFS 
junction is easily realized based on the current fabrication 
technology ji^^ 



VI. CONCLUSION 

We have developed a quantitative theory, which describes 
the properties of the DOS and the current-voltage character- 
istics of a SIFS junction in the dirty limit. We considered the 
case of a strong insulating barrier in a SIFS junction such that 
the left S layer and the right F/S bilayer are decoupled. In 
this case we can obtain the current-voltage characteristics of 
a SIFS junction in the framework of standard tunnel theory. 
In order to calculate quasiparticle current we first calculated 
the DOS in the ferromagnetic layer of the F/S bilayer. We 
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described the DOS behavior as a function of parameters char- 
acterizing properties of the ferromagnetic layer. In our theory 
we consider three such parameters: thickness of the ferromag- 
netic layer df, exchange filed h, and magnetic scattering t^. 
We have discussed the DOS properties paying special atten- 
tion to the DOS oscillations at the Fermi energy. We have 
also proposed to measure the exchange field in experiments 
on weak ferromagnets by measuring the DOS peak atE = h. 
We compared the results, obtained with a self-consistent nu- 
merical method, with a known analytical DOS approximation, 
which is valid when the ferromagnetic layer is thick enough. 

Using the numerically obtained DOS we have calculated 
the current- voltage characteristics of a SIFS junction and have 
observed features which are the signatures of the proximity 
effect in the S/F bilayer. We showed that there exists typical 
shape patterns of current- voltage characteristics related to the 
typical DOS structures in the ferromagnetic interlayer. 

Finally, we have calculated the macroscopic quantum tun- 
neling escape rate for the current biased SIFS junctions by 
taking into account the dissipative correction due to the quasi- 
particle tunnehng. Based on this we concluded that the influ- 
ence of the quasiparticle dissipation on the macroscopic quan- 



tum dynamics of SIFS junctions is small, which is a great ad- 
vantage of SIFS junctions for qubit applications compared to 
other types of ferromagnetic n - junctions. 
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